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Random bond Ising systems on a general hierarchical lattice are considered. The 
inequality between the specific heat exponent of the pure system, a p , and the crossover 

exponent 0, a p < 0, gives rise to the possibility of a negative a p along with a positive ^, 

leading to random criticality in disagreement with the Harris criterion. An explicit example 
where this really happens for an Ising system is presented and discussed. In addition to that, 
it is shown that in the presence of full long-range correlations the crossover exponent is 
larger than in the uncorrelated case. 



One of the most famous and important results in the study of systems affected by 
quenched disorder is the Harris criterion [1]. The widely accepted form of the criterion, is 
that in ferromagnetic systems with random interactions the randomness is irrelevant if a p , 

the specific heat exponent of the corresponding pure system, is negative, while for systems 
with positive a p the random system exhibits different critical behavior. Specifically, ^ = a p , 

where ^ is the crossover exponent from pure to random criticality. It is surprising, however, 
that in spite of the long time since its proposal and quite a number of alternative derivations 
[2-10], a rigorous proof of the Harris criterion is still lacking. In fact, the only rigorous 
results in this field [11,12] yield information about the exponents of the random system 
(cc r <0 or v r >2/d, d being the dimension of the system) without relating those to the 
exponents of the pure system. Some years ago, a counterexample to the Harris criterion was 
presented for Potts models on a certain hierarchical lattice (HL) [13]. It was shown that, for a 
g-state Potts model, it is possible to find a window of rather large q's in which a p is negative 

but the disorder is relevant. It is unclear whether the Harris criterion is also violated in the 
Ising model, because in contrast with high q Potts models [8] it is quite difficult to obtain 
negative a p 's in Ising systems. 

The above situation motivated the present work, which yields some rigorous results 
regarding the Harris criterion for the Ising model. Because on regular lattices any attempt to 
obtain some rigorous results via a renormalization group (RG) procedure at the pure fixed 
point is hampered by the generation of correlations, I chose to work on general HLs. This 
produces exact results because correlations are not generated in such systems under RG (note 
that some derivations of the Harris criterion [8-10] assume the same in one way or another). 
The results I obtain here are the following: 

(1) 1 define two random interactions to be fully correlated if both are identical in each 
representation of the randomness. I define the full correlations as long ranged if they survive 
the renormalization procedure for any number of steps. In the presence of full long-range 
correlations, the crossover exponent, to be denoted as ^ corr , is larger than 0, the crossover 
exponent of the uncorrelated system. 

(2) Explicit examples are constructed where a p is negative but still the randomness is 
relevant. This is true for uncorrelated and the correlated systems. 
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(a) (b) 

Fig. 1. A general HL is described schematically. In (a), the shaded area consists of a set of lattice 
points, where some of the pairs are joined by bonds GC,/3,.... In (b), a typical shaded area is 
represented. The full lines are bonds to be iterated in constructing the lattice. 

Consider a general HL described schematically in Fig. 1. The shaded area shown in (a) 
consists of a set of lattice points where some of the pairs are joined. In (b), atypical shaded 
area is represented. The lines are bonds to be iterated in constructing the lattice. All bonds 
carry a coupling J 12 , governed by a distribution P( J a ) that is identical for all bonds. The 
renormalized coupling is given by 

Ju=fV*}, (1) 

where / depends only on couplings associated with the pair of sites (1,2) (the shaded area, 
Fig. 1). This implies that J tj and J lm are not correlated if the pairs (ij) and (l,m) are not 

identical. The renormalized distribution P(J) is given by 

Hl 2 ) = \Y{dJ l a 2 P(J l a 2 )S[J l2 -f{J 12 }] (2) 

a 

and may serve to derive an infinite set of equations for the renormalized moments. Let us 
denote 

*]={{&„)'), (3) 
where SJ a denotes the departure from the pure fixed point J l a 2 = J* for all a . The recursion 
equations for the moments read 

/>G,[/],r 2 ,...]. (4) 

The pure ferromagnetic fixed point is assumed at J* > 0, so that 

r'=0. (5) 
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In a recent paper [14], the following results were proved by considering the matrix 

(drJdr^OA...). 

(a) The eigenvalues of the matrix, X { , are given by 

4=E(/ a )'' where f a =JL ( r_ r) . (6) 

The sum is over all n bonds a of the rescaling volume associated with the pair (1,2) (the 
superscript "1,2" is omitted hereafter) and the partial derivative is taken at the point where all 
those couplings equal J* . 

(b) All the eigenvalues are positive. 

(<M- + i<4- (7) 

As a consequence of (c), the leading crossover exponent is 

^- ln ^ (8) 



z In 4 
and 

^« P , (9) 
where a , the specific heat exponent of the pure system, may be easily obtained by 
consideration of the free energy per bond, 

2-a p = — . (10) 
ln/ij 

The equality sign in (9) holds only if all the bonds in the shaded area of Fig. 1(a) are 
equivalent. 

We turn now to the case where full long-range correlations are present in the system. 
In the following, I define full long-range correlations for a general HL. Consider a HL, 
generated using some specific generator such as the one presented in Fig. 2(a), where as 
before, the different J a 's (a = \,...,5 here) represent the independent random interactions 

attached to the different bonds. In Fig. 2(b), some of these bonds are taken to be fully 
correlated, namely, bonds a and j5 are fully correlated if, say, J ' p = J a so that indeed there 

is only one independent random variable, say J a . In order for this property to survive the 

renormalization procedure, the bonds have to be distributed in such a way that this property 
is true on all scales [as demonstrated in Fig. 2(c)]. In this way, correlations are nonlocal at all 
scales and thus long ranged. 
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(a) 



(b) 



(c) 



Fig. 2. In (a), an example of a HL, generated using some specific generator, is given. The different 

J a 's (OC = 1,...,5 here) represent the different random interactions attached to the bonds. In (b), 

some of the bonds associated with that generator are taken to be fully correlated, that is, the 

interactions J a attached to these bonds are taken to be identical. This must hold at all scales in 

order for the correlations to be nonlocal and survive the renormalization procedure. In (c), we 
zoom in on one renormalization step of a possible configuration that generates such correlations. 
The bonds are distributed in such a way that the initial setup of correlated bonds defined on the 
generator shown in (b) remains true on all scales. In this way, correlations are nonlocal on all 
scales and thus long ranged. 

Let us assume then that the n bonds inside a rescaling volume are divided into m 
subsets Q t , with t = l,...,m and \<m<n , each containing n t fully correlated bonds, so that 




n t =n. We can write 



because all the couplings in Q t equal J t . Finally, we use Eq. (6) for the X { 's to write 



This is true since f a >0 for all a 's [14]. The inequality <j) < $, orr thus immediately follows 

and the equality sign holds only if no correlations are present. A special case of the above 
inequality was obtained a long time ago by Andelman and Aharony [10]. Using the Migdal- 
Kadanoff (MK) renormalization scheme, they considered d- dimensional cubic systems with 
quenched bond disorder that are correlated along d 1 dimensions and showed that the 
crossover exponent and thus the Harris criterion are modified, respectively, to 
corr = cc p +d l v p and 2- (d -d l )v p < 0. Since the MK transformation performed on a 

regular lattice corresponds to an exact RG transformation performed on a diamond HL 



fUJ = gU t }, 



(11) 



m ml \ n 



^ OIT =Zfe)'=Z I A >£(/*)'' =4 (uncorrelated). 



(12) 



t=\ t=\ \a<Ef2 t J a=\ 
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[15,16], for which all bonds are equivalent and thus a = </> , their result reads 



fcort = <f> + d,V > . 

In the following I present an Ising system, for which a p is negative while <j> 2 is 

positive. Consider the HL generated using the generator presented in Fig 3. It is made of p 
parallel branches, each of which contains a single site, connected on one side with a single 
bond and on its other side with q parallel bonds. The rescaling volume is then n = p(q + 1) 

and, denoting fU a = K a , the renormalized coupling is given by 



K = f{K a ,K ap ) = \±\n 



cosh 



cosh 



(13) 



where K* is determined by 

K* 



\f ln cosh^ + DiT] 



(14) 



Taking the first partial derivatives of Eq. (13) with respect to K a and K a/3 , and using Eq. 
(6), the following expression is obtained for the eigenvalues; 

+ l)K*\- toW, - » K1 ) < + fa*, + D^-] + «anh [( , - I)K"]) ■}. (15) 

a p and 2 are next calculated according to Eqs. (8) and (10). Numerical results for K* , a p 
and 2 are given in Table 1 for different values of p and q. Note the cells emphasized in 
gray, where the Harris criterion mistakenly indicates pure criticality. 





(a) 



(b) 



Fig. 3. A HL generator is shown that is capable of demonstrating examples for the Ising model for 
which the Harris criterion fails to predict random criticality. It is made of p parallel branches, each 
of which contains a single site, connected on one side with a single bond and on its other side with 
q parallel bonds. In (a), the p parallel branches are shown. Each of the bold bonds represents a set 
of q parallel bonds as explicitly drawn in (b). 
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I also considered an example for the case where some of the bonds are fully correlated. 
For simplicity I have chosen the case in which the q parallel bonds of each of the p branches 
are fully correlated. Therefore, in each branch, we may replace the q bonds of Fig. 3(b), each 
carrying a coupling J aj3 , with an effective coupling qj' a , so that now 

\oMK a +qK' a ) 



cosh(K a -qK'J 



(16) 



and 



X- = f {(tanh[(? + \)K*] - tanhfa - l)iT ]) '- 



(17) 



+ q' (tanh[(? + \)K* ] + tanh[(? - \)K*]) ' } . 
While K* , \, and therefore a , are left unchanged, and therefore ^ corr are now 
increased. Numerical results for con . are also given in Table 1 . 



p 


2 


3 


4 


q 


1 


2 


3 


1 


2 


3 


4 


1 


2 


3 


4 


5 


K* 


0.609 


0.281 


0.185 


.361 


.175 


.116 


.087 


.261 


.128 


.085 


.064 


.051 






























-.676 


-.111 


.224 


-.902 


-.478 


-.090 


.117 


-1.19 


-.762 


-.319 


-.077 


.069 






























-.676 




.339 


-.902 


-.379 




.203 


-1.19 


-.670 


-.220 


.006 


.138 


Aon 


-.676 


.778 


1.26 


-.902 


.345 


.919 


1.16 


-1.19 


.040 


.680 


.961 


1.11 



Table 1. Numerical results for K , a <f> , and ^ COIT for different values of p and q with respect 
to the HL shown in Fig. 3. Note the cells emphasized in gray, where CC < < (j) so that the 
Harris criterion mistakenly indicates pure criticality. Also note that always CC < (j) < (j) cors . 



To conclude, I wish to discuss the possible relevance of my results to regular lattices. 
One should note that the two properties discussed here, that of bonds inside a rescaling 
volume not being all equivalent and that of full long-range correlations being present, each 
by itself makes <j> 2 > a p possible. This may suggest, at first sight, that starting on a regular 

lattice with a set of nonequivalent interactions, say, nearest neighbors and next-nearest 
neighbors, may violate the Harris criterion. On HLs we see, however, that it is also important 



7 



that these properties are invariant under the RG. On regular lattices, though, we expect that 
the detailed structure of nonequivalent bonds will disappear under RG and therefore the 
Harris criterion is expected to hold. Cases where on a regular lattice we may find negative 
a p but still randomness is relevant must thus be such that the nonequivalence of bonds does 

not disappear under RG. This can happen in principle when the interactions are algebraically 
long ranged or when the correlations of short range interactions are long ranged. 
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